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Abstract 

The invariant mass distribution of the di-photons from the decay of the lighter scalar Higgs boson(h) to 
be carefully measured by dedicated h search experiments at the LHC may be distorted by the di-photons 
associated with the squark-gluino events with much larger cross sections in Gauge Mediated Supersymmetry 
Breaking (GMSB) models. This distortion if observed by the experiments at the Large Hadron Collider at 
7 TeV or 14 TeV, would disfavour not only the standard model but various two Higgs doublet models with 
comparable h - masses and couplings but without a sector consisting of new heavy particles decaying into 
photons. The minimal GMSB (mGMSB) model constrained by the mass bound on h from LEP and that on 
the lightest neutralino from the Tevatron, produce negligible effects. But in the currently popular general 
GMSB(GGMSB) models the tail of the above distribution may show statistically significant excess of events 
even in the early stages of the LHC experiments with integrated luminosity insufficient for the discovery 
of h. We illustrate the above points by introducing several benchmark points in various GMSB models - 
minimal as well as non-minimal. The same conclusion follows from a detailed parameter scan in a simplified 
GGMSB model recently employed by the CMS collaboration to interpret their searches in the di-photon -|- 
/7t channel. Other observables like the effective mass distribution of the di-photon -I- X events may also 
reveal the presence of new heavy particles beyond the Higgs sector. The contamination of the h mass peak 
and simple remedies are also discussed. 
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1 Introduction 



The search for the Higgs boson/bosons and the study of their properties tops the priority Ust of the on going 
experiments at the LHC at 7 TeV as well as of the upcoming experiments at the highest attainable energy of 14 
TeV. This will shed light on the mechanism of electroweak symmetry breaking (EWSB) and the generation of 
the particle masses. The Higgs sector of the very successful standard model (SM) of particle physics responsible 
for EWSB is indeed very simple. It consists of a single neutral, scalar Higgs particle. However, this economical 
Higgs sector turns out to be the Achilles' heel of the SM. 

The scale of EWSB is expected to be around 100 GeV. On the other hand the naturalness or the hierarchy 
problem [1] [5] comes into play, should there be some new physics beyond the SM characterized by a much higher 
energy scale. The very existence of the new scale tends to push up the EWSB scale far beyond the expected 
magnitude. For example, the Higgs boson mass blows up to values closer to the new scale, in contrast to the 
preferred range [3] ruh = 89^26 suggested by EW precision data (the hierarchy problem), unless the parameters 
of the SM are extremely fine tuned (the naturalness problem). Even the absence of any compelling evidence in 
favour of a new physics model does not allow us to overlook the above issues, since one cannot wish away the 
Planck scale, the scale of gravity which is always there. 

The above theoretical inconsistency has led to several extensions of the SM. The minimal supersymmetric 
standard model(MSSM) pP removes the above problem elegantly in the limit of exact supcrsymmetry (SUSY) 
and stabilizes the Higgs mass. However, SUSY must be broken by the soft breaking terms, since the negative 
results of sparticle searches at various colliders indicate that the superpartners (or the sparticles) must be 
considerably heavier than the corresponding particles. In spite of SUSY breaking, the naturalness problem 
remain under control, provided the masses of the sparticles are ^ 1 TeV [4 or smaller. 

In the MSSM the Higgs sector must necessarily be extended and the spectrum consists of two neutral scalars 
{h and H), one neutral pseudo scalar {A) and two charged scalars {H^) [i, 5 . Enormous amount of work has 
been done in developing the strategies for searching these bosons at the LHC [51 [7] . Adding any new dimension 
to the above studies is not the aim of this paper. 

However, at the LHC these Higgs bosons are likely be produced along with a variety of superpartners of 
the standard model particles - the sparticles. In particular the strongly interacting sparticles may indeed be 
produced in numbers much larger than that of typical Higgs induced events. Even if a small fraction of such 
events pass through the selection criteria for the dedicated search of any of the above Higgs bosons, then the 
shape of certain distributions expected for Higgs production alone may change significantly, indicating that 
the Higgs boson may be part of a framework larger than the SM. Such an unexpected shape is likely to be 
prominent near the tails of the distributions under study where contributions from the Higgs signal as well as 
the SM background are naturally small. In this paper we shall illustrate this possibility with a specific example. 
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Some of the above issues have recently been addressed in the context of charged Higgs search by a counting 
experiment during the LHC-14 TeV run |5]. It has been demonstrated that both in the unconstrained MSSM 
and in the minimal supergravity(mSUGRA) model [5], SUSY events can seriously affect the size of the charged 
Higgs signal unless additional selection criteria are introduced. On the other hand, the combined squark-gluino 
and charged Higgs events may establish physics beyond the standard model (BSM) physics at a confidence level 
higher than that attainable by the Higgs signal alone. This enhancement would disfavour not only the SM but 
also its extensions with two Higgs doublets. After this important hint regarding BSM physics, a clean charged 
Higgs sample may be recovered by additional kinematical cuts. It needs to be emphasized that the above issues 
are relevant in any extension of the SM with a new particle sector as well as an extended Higgs sector. 

In this paper we focus our attention on the search for the lighter CP even scalar Higgs boson(h). The h boson 
in the MSSM has properties very similar to the SM Higgs boson in a large volume of the parameter space (the 
so called decoupling region). Thus the discovery of a single neutral Higgs scalar, which is a likely scenario 
in the early stages of the LHC experiment, will yield very little information regarding the underlying theory. 
The situation may change dramatically if SUSY is realized in nature in such a way that a sizable number of 
sparticle induced events change the characteristics of the expected h signal significantly. This can convincingly 
establish the BSM origin of the discovered boson. Even the discovery of one or more heavier Higgs boson with 
no counterpart in the SM, would indicate only an extended Higgs sector but not the existence of new heavy 
sparticles. 

The mass of h is always bounded from above (mh < 140 GeV) 4j in supersymmetric theories. As a consequence 
the search strategy for it is very similar to that for the standard Higgs boson with mass within the above bound 
especially in the decoupling region. If mh < 130 GeV, an important discovery channel would be the inclusive 
di-photon + X channel with a relatively modest integrated luminosity (£) of a few tens of fb^^ at the 14 TeV 
run of the LHC (see Figure 10.38 of [6 ). Here X stands for any other particle accompanying the h signal. This 
signal stems from h production at the LHC via several channels, followed by the one loop decay h — ?> 77. The 
suppression of the signal due to the tiny branching ratio (BR) of this decay ( ^ 10~^) is adequately compensated 
by the relatively low SM background. In the MSSM also a small peak in the 7 7 invariant mass distribution 
around the h mass would establish the h signal (see Figure 11.37 of [B]). 

Among various models of softly broken supersymmetry the Gauge Mediated Supersymmetry Breaking (GMSB) 
models (TUl |TT] (to be briefly reviewed below) predict a large number of di-photon + X events at LHC energies 
coming from the production of strongly interacting sparticles. In this paper we wish to focus on the impact of 
these events on the di-photon invariant mass distribution to be studied by the dedicated h-search experiments 
with utmost care. 

In these models the gravitino with negligible mass turns out to be the lightest supersymmetric particle (LSP) 
and, hence, the carrier of missing energy. If the lightest neutralino {xD is the next lightest supersymmetric 
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particle (NLSP), it can decay into the LSP and a photon with a large BR. In R-parity conserving supersymmetry, 
sparticles are produced in pairs and each of them eventually decay, through a cascade, into a gravitino-photon 
pair. Each SUSY event in this scenario will, therefore, have the 27 + X topology. Thus the pair production of 
squarks and gluinos in any combination having a large cross section may leave their signature in the observed 
7 7 invariant mass distribution. 

The Plan of the paper is as follows. In the next section we briefly describe different GMSB scenarios and 
the present collider constraints on them. We also define our benchmark points for the analysis. In section 3, 
we present our main results on the di-photon invariant mass distribution in GMSB models at the LHC. Our 
emphasis will be on the ongoing LHC experiments. However, the possibilities for the 14 TeV experiments will 
also be briefly touched upon. Our conclusions will be summarized in section 4. 

2 Different GMSB scenarios, sparticle spectra and the current 
collider constraints 

To begin with we shall review the constraints from different collider experiments on the parameter space of the 
minimal GMSB (mGMSB) model [lOlin]. We also wish to add a few points regarding the approximations in 
obtaining these bounds and the resulting uncertainties which have not been sufhciently elucidated in the current 
literature. In this model the soft breaking parameters are generated via the gauge interactions of the SM. As 
a result the sparticles with the same SM quantum numbers but different flavours acquire the same mass. This 
keeps the potentially dangerous flavour changing neutral current induced processes under control. 

In the mGMSB model there is an observable sector consisting of the MSSM fields. The supersymmetry breaking 
sector consists of a gauge singlet chiral superfield Y, identified with the goldstino superfield in the simplest 
version of the model. The scalar and auxiliary components of Y are assumed to develop vacuum expectation 
values (VEVs) denoted by S and F respectively. In the early versions of this model the dynamics of generating 
these VEVs were not specified. The model also has a messenger sector consisting of superfields having 
the gauge interaction of the SM. The messenger superfields interact with Y via a tree level interaction in the 
superpotential. This generates a supersymmetric mass of 0{M) of the messenger fields, usually referred to 
as the messenger scale. A SUSY breaking mass squared splittings of the order F is also generated within the 
messenger supermultiplets. This SUSY breaking is communicated to the observable sector by gauge interactions 
between the messenger and the observable fields via higher order processes. A restriction in the messenger sector 
comes from the requirement of the unification of the coupling constants of the standard model. If the messenger 
superfields form complete GUT multiplets (e.g., 5 or 5 of SU(5)), the value of the unification scale or the 
GUT scale {Mg) does not change. However, coupling constant unification constrains the number of messenger 
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superfields(N) to be < 5. 

The simple mGMSB model is characterized by five parameters [lUJ [TT] 

A, M, N, tan (3 and sign(/i) 

where A = F / S is the SUSY breaking scale in the observable sector, M is the messenger mass scale, N is the 
number of messenger multiplets belonging to the 5 + 5 representation of SU(5), tan /3 is the ratio of the vacuum 
expectation values of the neutral Higgs fields in the observable sector and fi is the Higgs mixing parameter in 
the superpotential with magnitude fixed by the radiative symmetry breaking condition. It bears recall that the 
soft breaking trilinear (A) and bilinear (B) parameters are generated by higher order processes and are small 
at the messenger scale M. The weak scale parameters are then obtained by the standard renormalization group 
(RG) evolutions. The sparticle spectrum at the weak scale can be computed by both SUSPECT version 2.41 
[12] and IS A JET [13]. The phenomenology of this model has been discussed by various authors [14] 

The current experimental lower bound on the Higgs boson mass from LEP nih > 114.4 GeV 15 poses the 
strongest constraint on the parameter space of the mGMSB model. Given the input parameters one can compute 
the Higgs mass including higher order corrections. In this paper irifi > 114.4 GeV will be henceforth referred 
to as the stronger bound on the computed h-mass. However, there is an estimated theoretical uncertainty of 
about 3 GeV on the computed h-mass due to yet unknown higher order effects [16]. In view of this a point in 
the parameter space with computed nih > 111. 4: GeV ( the weaker h-mass bound) may still be acceptable. 
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Fi gure 1'. Variation of the computed Higgs mass with the SUSY breaking scale A (in TeV) for different tan f3=5 (red line) and 
20 (blue line) in the mGMSB model. No significant enhancement of the Higgs mass in the region tanl3=20 to 50 has been noticed. 
The horizontal lines correspond to the Higgs mass 111. 4 GeV and 114-4 GeV respectively. The top quark mass is fixed at 173 
GeV and M=2 A. 
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In Figure[T]we present rrih computed as a function of A. We take mt=173 GeV, A/=2A, A^=:l. It follows that 
for tan (3 = 5 (tan /3 = 20) the stronger h-mass bound yields a lower bound Amin > 257 TeV (118 TeV). If, 
however, the theoretical uncertainties are taken into account and the weaker to^ bound is used instead, relaxed 
bounds Amin > 161 TeV (79 TeV) are obtained. It is interesting to note that for nih ~ 120, A sa 300 GeV which 
corresponds to gluino mass (rrig) and average squark mass (mg) of approximately 2 TeV and 3 TeV respectively. 
Thus should the Higgs mass bound continues to be pushed upwards, the squark- gluino search at the LHC, as 
predicted by mGMSB, will be quite challenging. 
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Figure 2: Allowed parameter space in the mGMSB model in the A (in TeV)- tanf} plane. The other parameters are as in Figure 
1. In the shaded (pink) region, -fi is the NLSP. The magenta and blue lines correspond to the Higgs mass 114-4 GeV and 111.4 
GeV respectively. The horizontal black line indicates the latest neutralino mass bound set by D0 collaboration. 

The D& and CDF collaborations have put strong constraints on NLSP (x?) mass. In order to obtain quantitative 
results they consider SPSS slope in the mGMSB model which has only one parameter A. The other parameters 
are fixed and given by M =2A, tan /3=15, A^=l and fj, > 0. The lifetime of the Xi is not fixed by these parameters 
and it is assumed to be sufficiently short. It means that the photon coming from the decay of is prompt. The 
gaugino pair production processes are expected to be dominant at the Tevatron. The decay of two Xi produce 
two photons and gravitinos that give rise to missing energy. 

The latest constraint is from D0 collaboration who have obtained a lower bound 124 TeV on the scale A with 
6.3 fb^^ of data at Run II [ITJ. We have computed the leading order (LO) cross sections for the dominant SUSY 
processes at Tevatron energies (mainly electroweak gaugino pair production) using PYTHIA |18j . The next to 
leading order(NLO) cross section and the K factor with renormalization scale{fiii)= factorization scale (Ai_F)= 
2m~± is computed by PROSPINO [19l[20. Our cross section for A=124 TeV (assuming SPSS point) is in good 
agreement with the result in [17 . We find that the above bound corresponds to an upper bound on the SUSY 
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production cross section of 4.5 fb. For promptly decaying Xi this Amin yields m^o > 174 GeV. This bound 
is indicated in the A - tan/3 plane (Figure [5]) by the horizontal line. For comparison the stronger and weaker 
h-mass bounds are also indicated in the same figure. It follows that the DO bound on A,„i„ supersedes that from 
the stronger (relaxed) rn/i bound for tab /3 ~ 15 (6.5). It follows that irrespective of the uncertainties in rrih, 
due to higher order corrections and the choice of tan (3, values of A smaller than 124 TeV are disfavoured. This 
bound also supersedes the earlier CDF lower bound on SUSY scale 107 TeV|21|. The assumptions underlying 
this bound will be critically examined below. 

Strictly speaking the D0 bound is valid for the snowmass slope SPSS with variable A. The other parameters 
are M = 2 A, = 1, tan j3 ~ 15 and sign (/i) > 0. For N ^ 2, the scenario reduces to a f NLSP scenario 
which is not under consideration in the present paper. 

However, keeping other parameters fixed and varying tan /? we have checked that neither the cross sections of 
the dominant processes ( xt ~ X2 ^^"^ Xi ~ Xi P^-h production |17[|21| ) at the Tevatron energy nor the x~^ and 
X2 masses vary appreciably with tan /3. This is, however, expected since over the scanned parameter space the 
above gauginos are wino like to a very good approximation and have masses controlled by the SU{2) Gaugino 
mass M2 at the weak scale. The cross section, therefore, shows only a mild tan (3 dependence. 

The bino-like Xi mass on the other hand is determined by the U{1) gaugino mass parameter (Mi) at the 
weak scale related to M2 by the gaugino mass unification condition. We shall, therefore, use the D0 bound 
throughout this analysis unless the electroweak gauginos happen to be significantly mixed, which is often the 
case for non-minimal models to be discussed below. 
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Figure 3: Minimum allowed value of the SUSY scale A as a function of the top mass obtained from the stronger and weaker 
Higgs mass bound for two values of tanp . 

The neutralino mass bound from Tevatron has an important advantage. The predicted mh depends sensitively 
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Masses 


CASE I 


CASE II 


CASE III 




tan/3=20 


tan/3=20 


tan/3=20 




A=120 


A=120 


A=120 




M,„=2A 


Af™=5A 


M„ = 10A 


m~ 


997.6 


963.8 


959.1 




1320.4 


1293.6 


1280.3 




1264.4 


1235.8 


1221.2 




417.4 


420.3 


423.6 


"lei 


424.6 


427.5 


430.8 




209.7 


213.3 


216.4 




199.9 


202.2 


204.6 


m^o 


169.1 


162.8 


161.9 




321.4 


311.8 


311.1 


Xi 


321.1 


311.7 


310.9 




114.5 


114.3 


114.1 



Table 1: Variation of sparticle masses as a function of the Messenger scale (Mm) for fixed SUSY scale and tan /3 in the minimal 
GMSB model. 

on mt (recall the dependence in the radiative corrections to rn^)- The minimum value of A consistent with 
the weaker or stronger h-mass bound as a function of mt is presented in Figure |3] for two values of tan /?: 5 and 
20. The top masses on the X-axis are allowed by the current bound obtained by a recent CDF measurement 
mt = 173.13 ± 1.2 [ig. 

These curves illustrate the uncertainty in A„ji„ due to the present error in mt- The neutralino mass bound, on 
the other hand, shows no strong rrii-dependence. Throughout this paper we shall use mt ~ 173 GeV. 

In summary, while choosing the representative points for calculating the SUSY contribution to the di-photon 
invariant mass spectrum in the GMSB model, the Amin from the weaker h-mass bound (neutralino mass bound 
) should be used for tan ~ 6.5 ( tan > 6.5). 

Another theoretical uncertainty comes from the choice of the messenger scale M. In Table [T] we present the 
variation of the sparticle spectrum with M. We find that the variation of M within moderate ranges does not 
affect the spectrum and, hence, the collider signatures drastically. 

It should be emphasized that the CDF and D0 experiments actually constrain the quantity (Jsusy x BR{xi — > 
7 G)^, where asusv is the total sparticle production cross section within the kinematic reach of the Tevatron, 
BR (xi — > 7 G) is the branching ratio underlying each single photon event. The cross section limit quoted 
above is, therefore, subject to the model dependent assumption BR {xi 7 G)^ « 1. 

The above BR depends on the neutralino mixing angles and may differ significantly from 1, if the lightest 
neutralino is indeed an admixture of electroweak gauginos and higgsinos. The widths of the different decay 
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A(TeV) 




105 


0.950 


125 


0.907 


145 


0.874 


165 


0.849 



Table 2: Branching fraction of Xi iG for ta,n f3=20 and M= 2A in the niGMSB model. 

channels of Xi are given, e.g., in Eqs 3.7 - 3.12 of [TT]. The relative width of the Xi decay into 7 and Z is given 
by (following the notation of [llj ) 

r(x? ^ Z G) = {kz/ r(x? ^7 G) X (1 - (m^/m^o)2)2 

As long as the Xi is a pure Bino, which is almost always the case in mGMSB, the ratio {kz/ k-^)^ w tan^Ow- 
Moreover for relatively small A, m^o « mz, and the Z-channel is strongly suppressed compared to the 7 
channel. However, as A and m^o increases, the relative importance of the Z-channel increases even in the 
mGMSB model, due to the factor (1 — (m^/m^jro)^)^. As a result the BR {xi — > 7 G) becomes significantly 
smaller than one. This is illustrated in Table[21 For A « 125 TeV, which is the current limit from the Tevatron, 
BR (xi — 7 G) w 0.91. This will relax the limit on m^o by a few GeV only. However, as the Tevatron 
experiments become sensitive to higher A, the value of BR (x? 7 G) should be carefully folded into the 
calculation before extracting the m^o limit. In the following we shall also consider more general GMSB models 
where the value of the BR is considerably smaller than unity and must be taken into account before extracting 
any limit from the Tevatron data. 



BP I 


BP II 


A=165 TeV, M„= 2 A 
tan 13=5, Ni = N2 = Ns^l 
ruh =111.6 GcV 


A=125 TeV Af„= 2 A 
tan (3=20, Ni = N2 = N3 = I 
nih =114.8 GcV 



Table 3: Benchmark points in the minimal GMSB model. 



In view of the above discussions we have chosen the points in Table [3] which are consistent with the constraints 
discussed above, for estimating the squark gluino contribution to the 7 — 7 invariant mass distribution from 
di-photon + X events in the mGMSB model. 

For point I with tan /3 = 5, the computed is consistent with the weaker bound and A = Amm as obtained 
from Figure[2] This point yields gluino mass {rrig) = 1329 GeV and average squark mass {rriq) = 1725 GeV. For 
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Point II A = Amin as obtained from the m^o bound, since tan /3 = 20. This leads to nig = 1035 GeV and nig 
— 1334 GeV. The squark-gluino contributions to the di-photon + X events in these scenarios will be computed 
in the next section. 

In the mGMSB model a single scale (A) controls both the sfermion and gaugino masses. As a result these 
masses are of the same order. Over the years several extensions of the niGMSB models have been proposed. We 
shall collectively call them non-minimal models. For phenomenological studies such extensions can be realized 
by introducing new parameters in addition to the set of five parameters mentioned above. 

One possibility is that the gaugino masses may be suppressed with respect to the sfermion masses [23;. This 
scenario can be studied by introducing a parameter R (< 1) such that Aq = -RA5, where Ag(A5) is the scale 
in the gaugino ( scalar ) sector at the messenger scale. 

More recently theoretically well-motivated extensions of the mGMSB model reflecting the above feature have 
been constructed. In a wide class of models the gaugino mass scale (Ac) happens to be severely suppressed 
in comparison to the sfermion mass scale (A5)[24 . This suppression is a consequence of expanding around 
the lowest classical vacuum of the low energy effective theory [3S]. This can be avoided, e.g., if the vacuum is 
an excited metastable state. It is, therefore, possible to construct gauge mediation models where the ratio of 
gaugino and sfermion masses continuously vary from very small to large values |26j . 

The sparticle spectra and BRs for different R can be computed by using, for example, online version of ISAJET. 
For R < 1 the gluino mass is expected to be significantly smaller than that in the mGMSB model (R =1) with 
A = As- However, the chargino and the neutralino masses will also be correspondingly suppressed in the above 
model. Hence, the Tevatron bound does not allow arbitrarily small R as we shall see below. 

For example, with As=:124, which corresponds to the Tevatron lower bound in the mGMSB model (see Point II) 
, we obtain for R — 0.6, m^o = 101, m^o = 205 ~ "^^±; ni~—667 and BR {xi 'j G) =0.98. It is obvious that 
with such light electroweak gauginos, the combind electroweak gaugino production cross section at Tevatron is 
enhanced violating the bound. For tan /? = 5 (20) the modified bound turns out to be A5 > 203(202) TeV. The 
Higgs masses in this two case are 112.7(117.5) GeV which is above the weaker rrih bound. Thus in contrast to 
the mGMSB model, the lower bound on A comes from the Tevatron data for both low and high tan /3. It is now 
easy to check the lower bounds on the gluino and average squark masses. They are m~ w 1030 GeV and niq w 
2 TeV. Thus although the minimum gluino mass is less than or comparable to the corresponding values in the 
mGMSB model, the average squark mass is much larger. Thus in general the size of the total squark-gluino 
signal and, consequently, the SUSY induced di-photon signal is expected to be smaller in non-minimal models 
with R < 1 compared to the corresponding mGMSB model. 

It has further been pointed out that not only the overall mass scales in the gaugino and sfermion sectors may 
be different, the SU{3), SU{2) and U{1) contributions of messenger fields to the masses of the gauginos and the 
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388.55 
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355.00 


327.04 
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234.25 




788.82 


835.52 


724.92 


689.48 


688.44 


503.11 




-174.37 


-185.26 


-236.53 


-221.43 


-332.55 


-216.57 




-332.91 


-325.56 


-377.98 


-309.82 


372.92 


249.46 


m~o 


441.64 


375.01 


406.32 


323.65 


-398.81 


-286.55 


m~o 
X4 


-477.58 


-441.27 


-524.44 


-490.80 


-744.13 


-549.85 


m~± 

Xi 


-325.08 


-320.00 


-369.82 


-304.86 


-366.79 


-240.69 


m~± 

X2 


-481.08 


-443.14 


-524.45 


-484.62 


-735.66 


-543.03 


■mh 


111.38 


116.65 


111.01 


115.70 


111.02 


114.06 


rriH 


908.35 


860.75 


832.37 


711.81 


778.70 


516.90 


rriA 


907.39 


860.54 


831.35 


711.61 


777.63 


516.72 




911.02 


864.64 


835.29 


716.49 


781.82 


523.31 



Table 4: Spectra for benchmark points in different non-minimal GMSB scenarios. 

sfermions having different gauge quantum numbers, may also differ from each other. This happens, e.g, when 
the messenger fields do not belong to a complete multiplct of a GUT group (say, SU{b)) |23) . 

Let us concentrate on a phenomenological model characterized by 5 parameters in addition to M, tan /3 and 
sign(/z): R (defined above), A5, A^i, N2 and N^. In the limit R = 1 and A^i= N2= N3 = N, this model reduces 
to mGMSB. The explicit mass formulae for the sfermions and the gauginos and further references can be found 
in Appendix A of [53]. 

The sparticle spectrum in this model can be computed by ISAJET. The Bino, Wino and gluino masses at the 
messenger scale are proportional to Ni,N2 and A^3 respectively. Clearly if Ni,N2 > N3, the mass difference 



10 



between the electroweak gauginos and gluinos at the electroweak scale obtained by renormahzation group 
evolution, will be reduced. Thus in this scenario the Tevatron bound becomes compatible with relatively 
light gluinos and, hence, one obtains a larger number of SUSY induced di-photon events. We illustrate some 
representative sparticle spectra with several benchmark points ( BP III - BP VIII) in Table S) 

For each point in Tabled the dominant chargino-neutralino pair production cross section ( NLO ) as a function 
of A is computed with PROSPINO assuming m~± « m^o . The K-factors for some representative values of m~± 
are given in Table [S] Comparing the computed value with the the cross section upper bound from Tevatron 
data, we obtain Amm in each case. Each A5 presented in Tabled is above the corresponding h-min- 





K- factor 


100 


1.39 


200 


1.29 


300 


1.15 


400 


0.99 



Table 5: The K factors of 

xt'Xji production at the Tevatron for different m-.± . We choose mass of the final state as the QCD 

scale. 

It may be noted from Table|4]that for most of the non-minimal models of this type, the gluino masses are signif- 
icantly smaller than one TeV. Another point of phenomenological interest is the composition of the electroweak 
gauginos. In the non-minimal model ft as determined by the electroweak symmetry breaking condition often 
turns out to be comparable or even smaller than Mi. As a result the higgsino component of the NLSP changes 
significantly. In Table |6] we present Mi, AI2 and /i in each model. The BR {xi ^7 G) is also shown. It follows 
that whenever the NLSP develops a significant higgsino component, channels like Xi'^Z G, Xi ^ h G open 
up and BR {xi ~^ 1 G) decreases. This has been carefully taken into account in extracting limits from the 
Tevatron data as well as in computing the size of the di-photon signal at the LHC in different scenarios. 

General gauge mediated symmetry breaking (GGMSB) models have been proposed recently [2Z]- These theo- 
retically well-motivated models incorporate the features of the above phenomenological models. A large number 
of authors have recently studied the phenomenology of GGMSB [5H]. In such a general framework, all gaugino 
mass parameters are independent of each other and squark masses and gluino mass are not correlated. Thus 
squark masses can be lighter in contrast to the models discussed above. To illustrate this point, we take two 
benchmark points where squark masses are free parameters. In BP IX, wc take the point BP III and set all 
squark masses to be 1 TeV. The BP X point is identical to BP V with = rrig. 

The total SUSY cross section at the LHC in each scenario is also presented in Table El The leading order and 
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Points 




Mg 


Ml 


M2 


M 


Total NLO cross sec 


Br(xi + 


Effective cr-sec 


BP I 


1725 


1329 


240 


445 


884 


221 fb ( 19 fb) 


0.85 


160 fb ( 14 fb) 


BP II 


1334 


1035 


181 


339 


502 


1.1 pb ( 84 fb) 


0.91 


910 fb ( 70 fb) 


BP III 


1618 


601 


181 


349 


432 


10.5 pb (583 fb) 


0.89 


8.3 pb (462 fb) 


BP IV 


1509 


635 


192 


371 


363 


8.7 pb (414 fb) 


0.87 


6.6 pb (313 fb) 


BP V 


1495 


769 


247 


474 


397 


3.0 pb (112 fb) 


0.81 


2.0 pb ( 74 fb) 


BP VI 


1427 


736 


235 


452 


313 


4.1 pb (207 fb) 


0.79 


2.6 pb (129 fb) 


BP VII 


1439 


1081 


376 


718 


364 


725 fb ( 48 fb) 


0.38 


105 fb (7 fb) 


BP VIII 


1072 


811 


273 


525 


240 


4.8 pb (319 fb) 


0.35 


590 fb (39 fb) 


BP IX 


1000 


601 


181 


349 


432 


16.8 pb (822 fb) 


0.89 


13.3 pb (651 fb) 


BP X 


770 


770 


247 


474 


397 


11.0 pb (640 fb) 


0.81 


7.2 pb (420 fb) 



Table 6: The average squark mass, gluino mass, total NLO SUSY production cross section at the LHC for y/s=H TeV (7 TeV) 
and BR(xi -f + G) for benchmark points ( BP I to BP X). The elctroewak gaugino mass parameters Mi, M2 and /i are also 
given . Note that in case of BP I, BP II and BP VII, the EW gaugino production cross sections are are comparable to the 
strong production cross sections. 

next to leading order (NLO) squark gluino cross sections are computed by PROSPINO The relatively 

modest electroweak production cross sections are computed by PYTHIA and multiplied by the appropriate K 
factor, as discussed above. In each case the QCD scale is chosen to be equal to the mass or the average mass of 
the sparticles in the final state and CTEQ 5L and CTEQ 5M parton density functions have been used for LO 
and NLO cross sections respectively. The cross sections agree well with ref 

It is found that the scenarios with relatively light gluinos ( BP III - VI) have significantly larger cross sections 
compared to the representative mGMSB scenarios ( BP I and BP II). However, for R = 1, ( BP VII and 
BP VIII) the NLSP has substantial higgsino component. As a result the relevant BRs are also small in the 
corresponding models and, consequently, the number of the SUSY di-photon events is rather small as we shall 
see in the next section. 

Recently the CMS collaboration has reported negative results for squark-gluino search [3^ search in the di- 
photon + channel in the context of a simplified GGMSB model 31 . In this model the m (the NLSP 

Xi 

mass), Mq and Mg at the weak scale are taken as variables whereas all all other sparticle masses are fixed at 1.5 
TeV. The impact of the squark-gluino induced di-photons in this case on the 7 7 invariant mass data collected 
by dedicated Higgs search experiments in the context of this model will also be reported in the next section. 
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3 Footprints of squark-gluino events in the di-photon invariant 
mass data collected by the Higgs search experiments at the LHC. 



From the sparticle spectra presented in the last section we find that in all cases the pseudoscalar Higgs boson 
mass {Ma) is much larger than Mz- Hence we are in the deep decoupling regime. As a result the production 
cross sections and the BRs of the h-boson are practically the same as the SM Higgs boson having the same 
mass. For h-bosons with masses a few GeV above the LEP lower bound, the di-photon+ X channel is the most 
promising signal. 

As already mentioned in the introduction the prospect of discovering the h-signal in the di-photon+ X channel 
at the LHC, is not the main concern of this paper. This has been dealt with in great depth by the LHC 
collaborations pi ^T^ . Our main task is to study the possibility of distortion in the tail of the 7 — 7 invariant 
mass distributions in different GMSB scenarios in typical Higgs search experiments. For this we simulate by 
PYTHIA the SM backgrounds and squark-gluino events in the models discussed in Section 2 using the cuts 
employed by Higgs search experiments. 

We begin by following the standard selection procedures for Higgs search in the di-photon + X channel by the 
CMS collaboration [6] at 14 TeV. We note that ATLAS collaboration [32] has used same cuts both for 7 and 
14 TeV analyses. Thus our approach seems to be reasonable. 

We require events with exactly two isolated photons with pseudo-rapidity \ri\ < 2.5 and px greater than 40 GeV 
and 35 GeV respectively. We put the following isolation conditions on photons. 

1. No charged particles with px larger than 1.5 GeV/c should be present inside a cone with AR < 0.3 around 
the photon candidate where AR — ^/ATf + A^. 

2. The total Et of electrons or photons with 0.06 < AR < 0.35 around the direction of the selected photon 
candidate, must be less than 6 GeV in the barrel and 3 GeV in the endcaps. 

3. The total transverse energies of hadrons within AR < 0.3 around the photon candidate must be less than 
6 GeV in the barrel and 5 GeV in the endcaps. 

The irreducible SM backgrounds are from i)q 9—7-77 and ii) g g 77 events. We generate these backgrounds 
by PYTHIA 6.4.21 subject to the above cuts. The LO cross sections of the backgrounds are computed by 
PYTHIA. The K factors of these two processes are given in [B] for 14 TeV LHC(1.5 and 1.2 respectively). The 
NLO cross sections for the SM backgrounds at 7 TeV are not available in the literature. As a reasonable guess 
we take the same K factor for i) and (ii) to be 1.5 and 1.2. 
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The additional cuts related to photon reconstruction employed by the CMS collaboration cannot be implemented 
in our analysis with the toy detector of PYTHIA. We assume that their absence will affect the SUSY events 
and the background similarly and our main conclusions will be by and large valid. We have also ignored the 
instrumental backgrounds. We admit that the main purpose of this analysis is to illustrate the possibility of 
distortion in the shape of an expected distribution due to new physics and not to present very accurate numerical 
results. 

The sparticle spectra in different non-minimal models (BP III - VIII in Table H]) can be readily generated by 
online ISA JET IJ. The spectra in mGMSB (BP I and II in Table E]) and GGMSB models (BP IX and X ) 
have been generated by SUSYHIT [33]. Finally the spectra are interfaced with PYTHIA for event generation. 
For squark-gluino events the NLO cross sections are directly computed by PROSPINO (see section 2). 



Processes 


M^^ > 200 


M^-y > 300 


M-y^f > 400 


M^-y > 500 


M^^i > 600 


Born (SM bg) 


707.6 


219.5 


91.3 


42.8 


20.6 


Box (SM bg) 


78.9 


12.8 


2.5 


0.6 


0.1 


BP I 


5.5 


3.0 


1.5 


0.7 


0.3 


BP II 


18.2 


7.7 


3.1 


1.3 


0.6 


BP III 


106.6 


44.8 


17.7 


6.9 


2.8 


BP IV 


75.7 


32.0 


12.7 


5.0 


2.1 


BP V 


23.8 


12.2 


5.6 


2.5 


1.1 


BP VI 


38.4 


17.7 


7.3 


3.0 


1.2 


BP VII 


0.2 


0.1 


0.1 


0.0 


0.0 


BP VIII 


13.4 


6.1 


2.6 


1.1 


0.5 


BP IX 


134.5 


56.2 


22.0 


8.6 


3.3 


BP X 


128.9 


71.5 


36.0 


17.7 


8.6 



Table 7: Number of events for C = 1 ih ^ in different invariant mass bins in benchmark scenarios presented in Table 6. We 
consider bins with M-y-y > 200 Ge V or more at LHC 7 TeV run. 

We next compute invariant mass distribution of the di-photons from h-decays, the squark-gluino events and the 
SM backgrounds. We show in Tabic □ the total number of events in all bins with M.^^ > 200, 300, 400, 500 
and 600 GeV for different benchmark points for C = Ifb^^. We have checked that the number of events from 
the Higgs signal in these bins are indeed negligible but sizable contributions come from the other two sources. 
It is interesting to note that the contribution from the squark-gluino events are statistically significant in some 
cases. For example, the total number of background events (B) in all bins with Mj^ > 300 for C = 2fb^^ is 
464. The corresponding number of squark-gluino induced events (S) for in BP X is 143. Hence S/\/B = 6.64 
in this case. Similarly for BP IX the contribution from squark- gluino production is above 5 cr fluctuation of 
the expected SM background. The number of events for BP HI and IV are also reasonably large albeit with 
smaller significance at this C. 
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0.6 0.8 1 1.2 1.4 1.6 1.8 2 



M~ (TeV) 

Figure 4: The parameter space in a simplified GGMSB model (the coloured region) where the distortion in the 
tail of the 7 — 7 invariant mass distribution is observable at the ongoing LHC experiments with C <5 fb^^ (see 
text for the details). The parameter space below the red line has already been ruled out by the CMS collaboration 

Of course such a small integrated luminosity is insufficient for Higgs discovery via the di-photon + X channel 
at the on going LHC experiments 35 . However, it is encouraging to note that even at this low C, BSM physics 
can show up in the distorted tail of the 7 — 7 invariant mass distribution. Moreover, if observed, this would 
certainly disfavour the competing SUSY breaking models like mSUGRA which do not have natural sources of 
such di-photons . This observation would also disfavour a two Higgs doublet extension of the SM model with 
comparable h- mass and couplings. Such a model may predict the Higgs peak at the right place but not the 
appreciable distortion of the tail, appreciably. 

It may also be noted that the contribution of the mGMSB is negligible for values of A constrained by the 
Tevatron data. In contrast, some of the non- minimal scenarios contribute significantly. Moreover, for R — 1^ 
the SUSY contributions even in non-minimal scenarios are negligible. Thus if a distortion of the tail is indeed 
observed it may indicate a non-minimal model with R < 1. 



Mh 
(GeV) 


99 

NLO (pb) 


Vector Boson Fusion 
NLO (pb) 


Wh, Zh, tih 
NLO (pb) 


Br(h -J> 77) 


115 


39.2 


4.7 


3.8 


0.00208 



Table 8: Higgs production cross-sections in the SM at the LHC for y/s=lJ^ TeV and the branching ratio to the di-photon final 
state (taken from CMS TDR). 
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We next consider the simplified GGMSB model studied by the CMS collaboration [30], briefly discussed at 
the end of section 3. We fix m^o at 50 GeV. The di-photons originating from the squark-gluino events are 
studied for different combinations of the common squark mass Mq and the gluino mass Mg. The results are 
presented in Figure!?]. The coloured area corresponds to the region where S/^/B as defined above is > 5 for an 
integrated luminosity determined by the colour codes given on the vertical line next to the figure. In computing 
the significance we have considered the events accumulated in all bins with M^^ > 300. The parameter space 
already excluded by the CMS collaboration from the £ = 36 fb^^ data is shown by the region below the 
red line. Thus hints of GGMSB scenarios can be obtained for a fairly large parameter space. 

We next repeat the above analysis for experiments at LHC - 14 TeV. The cross sections in Table |S] are from 
the CMS TDR [6], Table 2.1 for rrih = 115 GeV at a proton-proton center-of-mass energy of ^/s^U TeV. The 
NLO cross sections of the dominant SM backgrounds are also taken from [B]. 



The irreducible SM backgrounds discussed before and the di-photon events stemming from squark-gluino events 
in different GMSB models are generated by PYTHIA. 




Figure 5: The di-photon invariant mass distribution for the SM backgrounds, the Higgs signal (nih =115 GeV) 
and squark-gluino events for the GGMSB benchmark point BP IX at the LHC- 14 TeV run with C = 1 fb^^. 

We present in Figure [S] the distribution of Af-^^ from the decay of a Higgs with to;i=115 GeV (the tiny blue 
histogram), the SM backgrounds (the black histogram) and the squark-gluino events from BP IX (the red 
histogram)for C = 1 fb^^ at 14 TeV center-of-mass energy. The combined distribution of the SM backgrounds 
and the SUSY events is also shown ( the purple histogram). It is found that the SM background falls off rapidly 
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Processes 


M^^ > zUU 


M-yy > 4UU 


Al~,^f > 60U 


T\ 'f ^ 1 r\r\r\ 

Mjj > lUOU 


Born (biVl bgj 


1275 


185 


54 


8 


rsox (bM bgj 


302 


15 


1 





BP I 


54 


21 


7 


1 


BP II 


214 


60 


18 


2 


BP III 


1876 


434 


98 


7 


BP IV 


1519 


358 


82 


6 


BP V 


547 


177 


51 


5 


rir VI 


621 


173 


42 


4 


BP VII 


2 


1 


1 





BP VIII 


157 


45 


12 


1 


BP IX 


2580 


551 


112 


7 


BP X 


1977 


653 


199 


22 



Table 9; Number of events for £ = 1 fb ^ in the invariant mass bins starting from M-y-y > 200 GeV at the LHC 14 TeV run. 

for higher values of the invariant mass. Whereas, the SUSY contribution is fairly constant for a significant part 
of that region. The distribution of events in different bins is shown in Table |9l 

On the basis of the earlier simulations one does not expect a statistically significant Higgs signal at this tiny 
C. Yet the distortion of the distribution far away from the Higgs peak is already worth noting. For example, 
the number (S) of the squark-gluino induced events for M^j > 400 GeV turns out to be 551. The corresponding 
number (B) for the SM background is 200. Thus S / VB, as defined above is 39. This suggests that the excess 
of events coming from sparticle production, if observed near the tail of the distribution, can not be dismissed as 
mere statistical fiuctuation. The number of events with Mjj > 200 GeV or more for other benchmark points 
are given in Table IHl 

0.25 I 1 1 1 1 1 1 




3000 



M (GeV) 
eff 



Figure 6: Normalized effective mass distribution of the Higgs signal (Red dotted line), the GMSB benchmark 
point BP IX (blue) and the SM di-photon background (green). 



17 



The presence of squark gluino events may be revealed by other distributions as well. We present in Figure |6] the 
effective mass distribution of the Higgs signal , the SM background and the SUSY events for BP IX for £ = 
1 fb~^ at 14 TeV. We have defined the effective mass Me// as the scalar sum of jet pr's, lepton pr's, photon 
Pt's and missing transverse momentum: 

Meff ^^PT + ^ PT+ ^ PT+ pIt- (1) 

jet lepton photon 

The peak of the distribution at a high M^ff strongly indicate some BSM physics with new heavy particles and 
disfavour models with an extended Higgs sector but no new heavy particles. 




(TeV) 

Figure 7: Same as Figure^ hut for LHC experiments at 14 TeV with C <1 fb^^. 

In Fig [7] we present our results for the simplified model [30] . The conventions are the same as in Figure 2] except 
that we focus on the distortion of the M^^^ distribution for C < 1 fb^^ which will be accumulated at a very 
early stage of experiments at 14 TeV. A large region of the parameter space may indicate the distortion looked 
for. 

A few di-photon events resulting from sparticle production in some GGMSB model may enter into the Higgs 
peak. In Table [TU] we present the di-photon + X events in the neighbourhood of the Higgs peak at 115 for 
different scenarios for C = Ifb^^ at 14 TeV. In a bin of rrih ± 1.4, as suggested in we find the number of 
genuine Higgs induced events to be 43. The SUSY contributions to the same region are 33, 24, 47 and 19 for 
BP HI, BP IV, BP IX and BP X respectively. Thus the presence of SUSY events is indeed significant. 

Of course a suitable cut can easily eliminate these unwanted events. For example, the distribution in Figure [6] 
suggests that a cut of M^ff < 300 can eliminate the SUSY contributions around the Higgs mass peak, while the 
Higgs signal and the SM backgrounds remain practically unaffected. Before the cut, however, the non-standard 
origin of the Higgs boson will be revealed by the Mg// distribution. The SUSY events around the Higgs peak 
can also be eliminated by putting an upper cut on the ^t- However, as already discussed the M^ff distribution 
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Processes 


109 


110 


111 


112 


113 


114 


115 


116 


117 


Born (bM bgj 


90 


87 


83 


82 


81 


85 


76 


72 


67 


Box (biVl bg) 


66 


66 


60 


59 


57 


54 


57 


49 


50 


BP I 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


BP II 


1.2 


1.2 


1.1 


1.2 


1.2 


1.3 


1.2 


1.2 


1.1 


BP III 


10.6 


11.4 


11.9 


11.7 


11.1 


12.2 


11.8 


11.5 


11.3 


BP IV 


8.1 


8.1 


8.7 


8.1 


8.1 


8.8 


9.0 


8.4 


8.8 


BP V 


1.8 


1.9 


1.9 


2.1 


1.9 


2.0 


1.6 


2.0 


2.1 


VI 


2.1 


2.6 


2.7 


2.7 


2.8 


2.3 


2.7 


3.4 


2.5 


BP VII 


0.0 


0.0 


0.2 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


BP VIII 


0.6 


0.8 


0.8 


0.7 


0.8 


2.7 


0.5 


0.6 


0.9 


BP IX 


16.6 


19.6 


16.5 


16.8 


17.1 


17.6 


16.8 


17.1 


17.5 


BP X 


6.0 


6.5 


6.3 


6.1 


6.8 


6.7 


6.5 


6.5 


6.9 



Table 10: Number of events in the invariant mass bins 109 to 117 GeV from different sources after all cuts for C = 1 ih ^ at 
the 14 TeV run. The total Higgs cross section after passing all cuts is around 43 fb leading to 43 events at 1 ih~^ distributed in 
114- 116 GeV bins. 

provides a strong hint for the presence of new heavy particles, whereas large /Et can come from massless 
particles and/or jet energy mis-measurements. 

4 Conclusion 

In this paper we analyze the invariant mass distribution of the di-photons to be minutely studied by the 
dedicated Higgs search experiments in the di-photon + X channel at the LHC. In the SM only a tiny peak 
at rtih above the SM background is expected. In GMSB models, in contrast, the di-photons stemming from 
squark-gluino events may distort the tail of the distribution far away from the Higgs mass peak revealing the 
non-standard origin of the Higgs boson. In this paper we explore this possibility both at LHC 7 TeV and 14 
TeV experiments by considering benchmark points in different GMSB models - minimal as well as non-minimal 
(Table 5)[I21III1I231I21[151[M1[13. We also consider a simplified model recently employed by the CMS 
collaboration to interpret their negative results for SUSY search in the GGMSB models and illustrate the above 
point through a detailed parameter space scan. 

We first consider the minimal GMSB [TOj [11] model subject to the lower bounds on rrih from LEP and on the 
lightest neutralino (NLSP) mass (or equivalently on the scale A) from Tevatron (see Figures [T] and . The 
resulting lower bounds on the squark-gluino masses are so strong that di-photon events induced by sparticle 
production with tiny cross-sections are unlikely to affect the h-signal. This has been illustrated with the help 
of two benchmark points (BP I and II, Table [3]) in section 3. 
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The D0 and the CDF experiments actually constrain quantity asuSY x BR{xi 7 G)^, where asuSY is 
the total sparticle production cross section within the kinematic reach of the Tevatron, BR (x? — 7 G) is 
the branching ratio underlying each single photon event. The bounds on A or the NLSP mass in the mGMSB 
model are derived with the assumption that the above BR is one. 

If BR {xi 7 G)^ is significantly smaller than 1, the limit on asusY and, consequently, on A or the NLSP mass 
become weaker. This can happen in the mGMSB model (Table [2|) and, more importantly, in its generalizations 
(Table 6). While extracting A„ii„ from the Tevatron data in different models we have taken this possible 
reduction into account and have selected our benchmark points accordingly (see Table 01). In some of the 
models the total squark-gluino cross sections are significantly larger than the expected values in the mGMSB 
model and the relevant BRs are not too small either (see Table [6]). More recently it has been shown that 
the distinct features of these phenomenological models can be incorporated into theoretically well-motivated 
modelsIMlEilSnilSI]- 

Our main results are as follows: 

1. The di-photon + X events from sparticle production may change the shape of the tail of the 7 — 7 invariant 
mass distributions which will be analyzed with utmost care by the LHC experiments. This distortion was 
shown to be statistically significant in several benchmark scenarios of the non-minimal GMSB models for 
both LHC 7 TeV (Table [7]) and 14 TeV (Table [9]) experiments. In both cases the traces of new physics 
beyond the Higgs sector may show up at the early stages of the experiments with C insufficient for the 
discovery of the h boson in the di-photon + X channel. Similar conclusions follow by analyzing a simplified 
model [20] ( Figures Hand [3). 

2. Other observables like the effective mass distribution of the di-photon events may also contain hints of new 
physics consisting of heavy particles (Figure |6l). Unexpected shapes of this distribution may disfavour not 
only the SM but i) extensions of it with larger Higgs sectors but no new heavy particles, ii) models with 
other SUSY breaking mechanisms like mSUGRA which do not have natural sources of di-photon events 
from sparticle production and iii) the mGMSB model subject to the constraints from LEP and Tevatron. 

3. The squark-gluino induced 7 — 7 events may contaminate the peak of the invariant mass distribution at 
rrih (see Table [TU| . To improve the purity of the h-signal one may implement cuts on suitable kinematical 
variables like an upper cut on the effective mass (see Figure IHl or an upper cut on the missing energy 
which eliminate the sparticle induced events in the Higgs peak. 
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